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Abstract 

in 

\ A Mathematica Notebook is presented which ahows for the transfer or any kind 

of polynomial expression to Matlab. The output is formatted in such a way that 
■ Matlab routines such as "Root" can be readily implemented. Once the Notebook 

has been executed, only one copy-paste operation in necessary. 
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^ I 1 Introduction 

^ ! Transferring a symbolic Mathematica expression to Matlab is a recurring issue 

O ! one for anyone using these two programs. On the one hand, the Mathemat- 

ica software is perfect to work out analytical calculations and derive involved 
symbolic formulas. On the other hand, Matlab offers a quasi-infinite and very 
^ ■ flexible environment for graphical, numerical or statistical treatment. The per- 

spective of working in both environment to solve a given problem is thus very 
attractive. As long as the equations involved remain simple, an expression 
derived in Mathematica can be straightforward reproduced in Matlab. Prob- 
lems arise for long and intricate formulas which cannot be easily copy-pasted. 
The transfer of a polynomial poses an additional problem: if one is willing to 
take advantage of the Matlab polynomial routines, the very definition of the 
function needs to be adapted because Matlab defines a polynomial from its 
coefficients. 



The Mathematica Notebook presented in this paper offers a simple solution 
to this problem. Given any polynomial expression derive in Mathematica, the 
Notebook generates a text which can be directly copy-pasted into Matlab. 
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2 Notebook description 



Consider a polynomial of the form, 

n 

Fix) =J2(^ii'^l- ■ ■Oik)x\ (1) 



where the coefficients ai...„ are themselves polynomial expressions of some 
parameters ai...fc. While both Mathematica and Matlah can deal with this 
kind of expressions, Matlab requires a different format if functions such as 
roots are to be exploited. More specifically, the polynomial P defined by Eq. 
([1]) will be expressed in Matlah under the form of a ID matrix, 

p = [a„(ai . . .ttfc) a„„i(ai . . .Ofc) ... 00(01 ... a^)] 

The Notebook which is about to be described generates the format needed by 
Matlah from the Mathematica expression. The starting point is the definition 
of the polynomial, 

In[l]:=P =a„(ai . . . afe)x" + a„_i(ai . . . + • • • + ao(ai • • • "fc); 

Note that the polynomial does not need to be specifically formatted like in 
Eq. ([1]) . The Notebook can perfectly deal with an expression where the coeffi- 
cients of the x*'s are not factorized. If the polynomial expression accounts for 
Greek characters (dealt with by Mathematica) , they will have to be replaced 
by some equivalents for proper Matlah treatment. For example, if the Math- 
ematica polynomial accounts for the parameters /3, 7 and u, the next line of 
the Notebook is, 

In[2]:=P Clear = P/. {/3 — )■ beta, 7 — gamma, u — )■ omega}; 

We then need the degree of the Polynomial, given by 
In[3]:='Deg = Exponent [PClear, x]; 

We now have Mathematica extracting the Oj's with, 

In[4j:=Coefs = Table[InputForm[FullSimplify[CoefRcient[PClear, x, 
k]]], {k, 0, Deg, 1}]; 

Here, the k^^ coefficient is first gathered through the Coefficient [PClear, x, 
k] command. The resulting expression is simplified with FullSimplify. Note 
that a simple Simplify can be inserted here if the argument is too involved 
for FullSimplify to execute successfully. 
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Cocfa . Table[IivutFann|TBlISuicrllt]'rCaeff±clBl[fe|7ClBK, x, n]]], {n, 0, Beg, 1)1; 
!>□ [Print , Deg + 2 - >i, ")=" , CDefs[[n]], " ; " ] , (n, I>eg + 1}] 

P(13)--(S2-^2* (Leta-"4'-glj-^Z* (gb * OnegaB -^4 - 3 ^ alpha-" 3 * DaiegaB-" 2 * S^-^ 2 + gb-"5--ZE--4 - 

2'alpha-gb-^2'-DmegaB-^2'- (2x' 2 + 2 * (1 + OiLegaB-^2) * 22 '2) + 2 alpha* gb '4^ 2s -"2 * (Zx-^2 t 2- (1 + OmegaB -^2) '■2s -"21 t 
2--gL-"3-Z2-"2* (2*alpba'2-- [2x'-2 + Zz'-2) + OmegaB-^2* (-1 + 2^alpha'2* (-15 t- ly-'-Z + 2e'2)))) t 
beta-" 2-- gb * ( (gb * OmegaB-" 2 - gb '3 * 2s -"2) '2 ^ ( (2 + OmegaE' 2) 'Zy.'-2 -i 2*(1 + CimegaB-^2) *2z' 2) + 

4--alpha'-2^QmegaB'-2^ (QmegaB'2-> {,2>.'-2 + 2*22-^2) - gb-"2*22-^2* ((1 + gb-"21 *2j4-^2 + 2*S2-^2)) + 
2 alpha - gb 3 ^ (Ix' 2 + 2 - Is -"2) ^ (- (OmegaB -^4* C-S -i- Sk-"2 + ) + gb -^2 * Is '2 (Ix -"2 + Zs-^2) + 

OmegaB'Ss- (2x-"2'' (-1 + gL'2--Ss-"2) + Zs-"2-- (-1 + gb'2^ (-15 + Zs'2n))) -t- 

+ .SS-"23 • H'aljiia' (OmegaB-^2 - gb-"2'Is-^2) * i:0megaB'2 - gb^£« (gb'-2-Sx-'2 + Sl-^^n -i- 

- fftMasE-t-S - (-1 + + Zs-"2) + gb '■ 4 * Is ■ 4 » (Zj( -" 2 + Zs-^ 2) + gb-^2 * Oaieg^-^: - Zs '2' (Ix -"2 - (-2 + gb'2 -Ss-"2) + Zs'-2^ 

t:-2 + *-^2* (-15 + + DiiegaE-^4' (Ix'^-- (1 - 2 - gb -"2* l2-"2) + Is'2-. tl - 2.ab^2 * (-6 + Zs'S))))))); 

.tBfi* [Zx"-^ *■ l^'Tj * {i4 4-alplie*Qsieg.aE-~4i2z'-£ - 4 alpha* gb'-4t DaegaB'-2 ( + (fi * Oma gaB 4 * (i 4- DaegaB'-^) Zz"2) * 

^•r aIliha*gb'-2»OnegaP'-2 * (OKf^q^E-^^' Zy-"2 - 2:fZE"4J - 4tgb'-3:>0aegaG-^2^S3''2»i;&:''2-k Zs-'a* Dmegffl-"2* (-l-.a('2*ZE-"2)) - 
gb'i-.EE^4T [a. f£e-'2 1- EE'-2) + (-10 ♦ 3^Zy.'-Z + 3:.Z!:-^2))J * 

Beta'4»- (^'-0B«S3B'4 - Ifi ^ alpha' 3 * OmegaB ■ 2 * 2z -" 2 + 3tqt-^5:>Es-'4» 4 » aipijai ^ 4 :> Zs 2^ (Zk '2 - 3 * (1 . OmegaB -" 2) :> Ee 2J - 
2=-alpha';gb'-2 *D»egaB-"2 :> (Zx-^2 + 4*(l-aBeijs^-"2)-.2l-'2) + 

4<.glj'-3^ZE'-2-. (alpha' 2^ (2* Zx' 2 + 3 .EE'-2i * Ctee-gqfi' 2 ^ [-1 ♦ alpha-"2 . (-10 + 2 ^ Zx' 2 + 3 ^ Zz -" 2) ) ) ) + 
Irsta'Si^ (^-"2, ((2 + QaegaB '2) ^Zx-"2 + 2»(1 ^ OasgaB-^S) - Zz '2) » ([)isegaB-"4 - 4:.^-"2*0sieg^-"2.Zz-^2 + 3^gh'4-.Zz'4) + 
4»^lI.t(a-"£- DiLeqaB-"2* (OmegaB-^2^ (Zx'2 + 2*Zs-"2) - 2-^-^2^Zs-^2^ ((1 + gb-^2) »Zx'2 + 2^Zs'"2)) + 2 * alpha ^ 3 * 
!£-gt^2.Zz-"2*(Zx'-4+4^Zx'2 + Zs-"2 + 3*Zz-"4) -□negaB'4^(ZK-^4+4*Zz-^2*(-l4-Zz-^2) +Zx'2*(-2 + 5*Zz'2)) *- DiiiegaB-^2 * 
/2k.'-4«' 1-1 + 2- iib-^2' Zs-^21 + 2 ■-Ss-"^ * (-2 + ah' 2-- f- 10 + 3- 2s '21 1 + Zx' 2 ^ 2s' 2 (-5 + 2 * ab-^2* (-5 + 4*2s' 21111 11 ; 

Fig. 1. Copy-pasting the Mathematica output to Matlah. 




The InputForm command is a key element of the transfer to Matlab. Mathe- 
matica can deal with symbolic fractions or power, but not Matlah. For example, 
the Mathematica expression. 



2l3 



c\t - u) 



(2) 



must be cast under the form 

a"2*b"3/(c"4*(t-u)), 

for Matlab to interpret it. The InputForm command is responsible for the 
Matlab formatting of the polynomial coefficients. Finally, the Table command 
ensures that the operation is implemented for all the a^'s, and the result is 
stored in the Mathematica array Coefs. 

We finally have Mathematica displaying the content of the array Coefs under 
a form which can be copy-pasted to Matlab, 

/n/^y.-=Do [Print ["P(", Deg + 2 - n, ")=", Coefs [[n]], ";"], {n, Deg + 

1}]; 



The command above thus prints all the coefficients of the polynomial from 
P(l), coefficient of x" to P(deg+1), coefficient of The output can then be 
copy-pasted to a ".m" Matlab file. Note that it is here necessary to "Copy As 
Plain Text" from Mathematica, as illustrated in Figure [1] 
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3 Conclusion - An example 



Two files have been attached to this article, as an illustration of the kind of 
transfer allowed by this technique. The Mathematica Notebook "Bridge to 
MatMab.nb" introduces a imposing polynomial which arises as the dispersion 
equation of a quite involved relativistic magnetized beam plasma system [Ij. It 
is a very involved polynomial in terms of the parameters x, Zx, Zz, a, jb, 13, fls- 
The degree of the polynomial in terms of x is 14. The expression results from 
the determinant of a dielectric tensor, and the coefficients of x are not gathered 
at all. The execution takes only a few minutes, which can be reduced even 
further by replacing the PullSimplify by the faster Simplify command in 
the coefficients calculation. 

After copying the output as shown on Figure [H the result has been directly 
pasted in the Matlab file "poly.m" and can be readily exploited in this Soft- 
ware. 
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